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We investigate the effect of piezoelectric (PZT) material on the flutter speed, vibra¬ 
tion mode and frequency, and energy harvesting power and efficiency of a flexible flag 
in various fluids. We develop a fully coupled ffuid-solid-electric model by combining 
the inviscid vortex sheet model with a linear electro-mechanical coupling model. A 
resistance only circuit and a resonant resistance-inductance (RL) circuit are consid¬ 
ered. For a purely resistive circuit, an increased electro-mechanical coupling factor 
results in an increased flutter speed, vibration frequency, averaged electric power and 
efficiency. A consistent optimal resistance is found that maximizes the flutter speed 
and the energy harvesting power. For a resonant RL circuit, by tuning the inductance 
to match the circuit frequency to the flag’s vibration frequency, the flutter speed can 
be greatly decreased, and a larger averaged power and efficiency are obtained. We 
also consider a model scale set-up with several commonly used commercial materials 
for operating in air and water. Typical ranges of dimensionless parameters are ob¬ 
tained for four types of material that span a wide range of solid density and rigidity 
values. We find that the resistance only circuit is more effective when the flag is 
placed in a lighter fluid (e.g. air), while the RL circuit is able to reduce the flutter 
speed when the flag is placed in a heavier fluid (e.g. water). 
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I. INTRODUCTION 


Piezoelectric (PZT) material has drawn enormous attention in recent decades due to its 
ability to convert mechanical kinetic energy into electrical potential energy, and vice versa. 
It has been applied to energy harvesting devices to power micro sensors/actuatorsii^, and 
to passively control structural vibrations^^-. The fundamental mechanism in the electro¬ 
mechanical process requires vibrations of the solid, which can be realized through natural 
flow-induced vibrations. 

A flexible elastic flag undergoes flow-induced vibration through the competition between 
the destabilizing effect of the fluid pressure and the stabilizing effect of the flag’s internal 
and bending rigidity. In particular, the flag will develop large-amplitude deformations if 
the flow speed is above a certain critical value, which we dehne as flutter. Much previous 
work has been conducted on this phenomenon through experiments^^—, inviscid and viscous 
numerical simulations!^^—, and linear stability analysis!^‘22i^. The critical flutter speed was 
found over a large range of mass ratio (ratio of solid to fluid inertial force), and increases 
as the mass ratio decreases. The effect of mass ratio on the flutter mode was observed 
in experiments and numerical simulations, and both found that the flag flutters in higher 
modes when the mass ratio is small and the hrst mode shape is never unstable for a cantilever 
platei^iSi. Post-flutter behavior was also discussed in previous work, and flags were found to 
transition from limit-cycle oscillations to more chaotic states as the fluid velocity increases 
beyong flutter— 

The objective of this work is to consider the effect of the piezoelectric material on the 
flutter speed, vibration mode and frequency, and energy harvesting power and efficiency 
of a ffexible ffag. Previous studies found that the electro-mechanical coupling effect of the 
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piezoelectric material acts to increase the bending stiffness of the system, and therefore will 
tend to stabilize the flag by increasing the critical flutter speecU^i^i. The output circuit also 
has an important impact on the flutter boundary. It is often described as a shunting circuit 
as the piezoelectric material is modeled as a current source in parallel with its internal 
capacitance^^i^^, and the output circuit provides extra paths for the current. A resistive 
shunting circuit with a single resistor and a resonance shunting circuit with resistor-inductor 
(RL) are two commonly used output circuits for energy harvesting and vibration control^. 
With the resistive shunting, Akcabay and Youngi2. studied the fluid-solid interaction of 
a bimorph piezoelectric cantilevered beam in a viscous flow using the immersed boundary 
method, and they found an optimal resistance value for energy harvesting efficiency. Michelin 
and Doare^^ studied the piezoelectric coupling effect on the local and global instabilities, 
as well as on energy conversion efficiency for a flexible plate using the double wake inviscid 
model. They found that the critical flutter speed increases as the piezoelectric coupling 
increases, and the energy conversion efficiency depends on piezoelectric coupling, loaded 
resistance, and the thickness ratio of the PZT patch to the substrate material. In another 
work, they studied bodies of small span using slender body theory^!. The increase of the 
flutter speed by the PZT material has an important impact on vibration control. On the 
other hand, higher effective system stiffness also indicates the flutter initiates at a larger flow 
speed, which will impact the energy harvesting potential and efficiency. In one aspect, this 
limits the effectiveness of the piezoelectric material in energy harvesting, especially in water 
and other heavier fluids because the required flutter speed is higher when the solid-to-fluid 
mass ratio is smaller. Some recent work suggests that this difficulty may be overcome by 
connecting the PZT flag to a well-tuned resonant shunting circuit. Xia et ai— found that 
the flag’s vibrating frequency can lock-in to the circuit resonant frequency for a parallel 
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connected RL circuit with a flapping PZT flag using the slender body theory. They showed 
a signiflcantly enhanced energy efficiency and a large decrease of the flutter speed in the 
lock-in regime. Similar results were obtained for the passive control of vibrating plates. 
Li et ai— investigated the effect of the RL circuit on the passive control of a hydrofoil. 
They also observed a maximum damping enhancement for an open-loop vibration control 
application via PZT with an RL circuit when the inductor is tuned to match the circuit 
resonant frequency with the foil resonant frequency. 

In this work, we extend the vortex-sheet model by Albert^ and include the piezoelectric 
effect by addition of a linear electro-mechanical coupling model to develop a fully-coupled 
model to predict the electro-fluid-solid response of a piezoelectric flag connected with an 
output electric circuit in incompressible and inviscid flow. In this model, the vortex layers 
induced by the body tend to sheets of infinitesimal thickness. A vortex-sheet model is 
efficient for slightly viscous flow because it can present the flow using vortex sheet, which 
is only one-dimensional in a 2D flow. It allows flow separations at prescribed locations 
(usually sharp edges), but it agrees quantitatively with Navier-Stokes models in predicting 
the strength of shed circulation, and the location of the primary vortices in the wake of an 
oscillated plated. Previous work also shows a good agreement of the vortex-sheet model 
in predicting the stability boundary of a flapping flagl^. Linear electro-mechanical coupling 
models were applied to cantilever flags for a single-patch PZT layer model^^i^ and a multiple- 
patch layer— model, and were found to compare well with experimental results at low to 
moderate strains. We then derive the important dimensionless parameters that govern the 
model, and investigate their effects on the critical flutter speed, vibration frequency, mode 
shapes and energy harvesting power and efficiency with the numerical simulation results. 
Next, we discuss the effectiveness and scalability of bimorph PZT flags for energy harvesting 
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and vibration control applications in air versus in water and for different substrate materials. 

The paper is organized as follows: Sec. II describes the vortex-sheet-plate model for the 
piezoelectric flag; Sec. Ill studies the effects of piezoelectric parameters on the stability- 
boundary; Sec. IV presents the scalability results for actual devices with different materials 
and fluids. The conclusions are presented in Sec. V. 

II. MODELLING 

A. Plate-Vortex-Sheet Model 

We consider here the motion of a cantilevered thin plate or beam in an inviscid flow. 
The beam is of chord length G, span length S, mass per unit chordwise length mg, and 
bending rigidity B. The flow and the motion of the beam is assumed to be 2D, so all the 
properties are uniform in the spanwise direction. The thickness of the beam is assumed to 
be much smaller than its chord and span, so we model it as a ID inextensible elastic sheet. 
Piezoelectric layers are combined in parallel on both sides of the beam by poling the layers in 
the same direction4^>2^. The beam moves under the pressure forces of a surrounding inviscid 
and incompressible flow of density (mass per unit volume) pf with an uniform upstream 
horizontal velocity U. A schematic figure of the piezoelectric beam is shown in figure [T](a). 

The instantaneous position of the beam is described as (C(s, t) = x{s, t) + iy{s, t), where s 
is the arclength. Assuming an Euler-Bernoulli model for the plate, the bending moment at a 
given position consists of an internal rigidity of the material and the piezoelectric coupling: 

M{s,t) = BK{s,t) — l3MV{s,t) (1) 

where B = El is the bending rigidity with E the material’s Young’s modulus and I its second 
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FIG. 1. (a). A cantilever beam of chord length C subject to a uniform flow of speed U and density p/. The 
arrows on the piezoelectric layer indicate the two layers are poled in the same direction, (b). Equivalent 
circuit parallel connected with a resistor, i? is the resistance, L is the inductance, Gp is the equivalent 
capacitance of the PZT material, I is the current generated by electro-mechanical conversion effect, 
moment of area, K{s,t) is the curvature of the beam, is a electro-mechanical coupling 
factor indicating the bending moment generated per unit voltage, and V{s,t) is the electric 
potential generated between the two piezoelectric layers. Therefore, the governing equation 
for ( is: 


msduC{s,t) = ds{T{s,t)s) - Bds{dsK{s,t)n) + {]Mds{dsV{s,t)n) - S[p]{s,t)n. ( 2 ) 


Here, T{s,t) is the beam’s tension force accounting for its inextensibility, and [p](s,f) is the 
pressure jump across the beam, s = dgC and h = is represent the unit tangent and normal 
vectors to the beam, respectively. The piezoelectric composite beam is parallel connected 
with a resistance-inductance (RL) output circuit. The voltage V{s,t) is governed by the 
linear dynamic equation for the circuit: 


dtV V 

CpdttV H-h — + = 0 

r I 


( 3 ) 


Cp 

where Cp = —, r = RC and I = LC are the surface density of the internal capacitance, 

(_y 
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resistance, and inductance, respectively with C as the chord length of the beam. The equa¬ 
tion is derived by representing the piezoelectric element as a current source in parallel with 
its internal capacitance2^>^ and has been used before in previous work by Xia et al.—. The 
equivalent circuit is shown in £gure[T]^b). Since the piezoelectric layer is parallel connected, 
the value of capacitance is twice that on each layer. The piezoelectric effect also plays a role 
in the trailing edge condition, where the pressure induced moment and forces are zero: 

T = = 0, /£ = (4) 

ID 

If we assume the voltage V = V(t) is constant along the piezoelectric layer, then the 
above equation can be rewritten by integrating along the arclength: 

CpdttV H — — —I" ^ j dttK.{s',t)ds' = 0 (5) 

In this formulation, /3Mds{dsV{s,t)n) = 0, equation ([2]) then reduces to the Bernoulli beam 
equation^^, i.e., without piezoelectric effect. In this case, the piezoelectric material affects 
the flag only through the boundary condition, and we can view the piezoelectric material as 
adding an extra concentrated moment at the free tip of the beam. 

The pressure forces are computed by the inviscid vortex-sheet model described by Albert^. 
This model allows flow separation at sharp edges only, but it agrees quantitatively with 
Navier-Stokes models in predicting the strength of shed circulation, and the location of the 
primary vortices in the wake of an oscillating plated. In the vortex-sheet model, the flow is 
computed in terms of the position and strength of a single vortex sheet. The vortex sheet 
consists of a “bound” part (Cb) on the beam, and it separates from the beam’s trailing edge 
into a “free” part (Cj). The position of the vortex sheet is denoted C('S,t), which is the 
same as the beam’s position when 0 < s < C. The strength of the vortex sheet is denoted 
7 (s, t), which is also the jump in tangential velocity of the beam. On the bound vortex sheet. 



the flow satisfles the no-penetration condition, which requires the body’s normal velocity to 
equal the flow’s normal velocity: 

Re(nd,C) = Re (n [u + j , c,M e C, (6) 

Here, the symbol denotes complex conjugate and U is the speed of the background 
flow. The special integral symbol in the equation denotes a principal-value integral, and 
the general solution 7 (s,t) has inverse-square-root singularities at the beam’s leading and 
trailing edge^. The last term is a regularized Biot-Savart integral using the smoothing 
parameter 5 of Krasny’s method^: 


where 


hs{s,t) 





ds', C(s, t) e Cf 


d(s,t) = do(l- exp (-jC(L,t) - C(s, f)|V45o)) 


( 7 ) 

( 8 ) 


with (5o = 0.2. This tapered regularization decreases the effect of the regularization term 
at the trailing edge as S{s,t) —?• 0 with s —)■ C, and it allows smooth vorticity and velocity 
dynamics for the free vortex sheet^. 

On the free sheet O/ in the wake, the vortex sheet evolves according to the Birkhoff-Rott 


equationr““: 


^ ^ ^ 2mJc,as,t)-as',t) 


—ds' + bs{s,t), C{s,t)eCf 


( 9 ) 


The pressure jump [p](s,f) across the beam is related to the vortex sheet strength 7 (s,f) 
using a version of the unsteady Bernoulli equatiort^^— : 


+ ds{{Tf{s,t) - T,{s,t))'j{s,t)) = — a,[p](s,f) (10) 

pf 

Here, Tf{s,t) and Ts{s,t) are the tangential velocities of the fluid and the solid respectively. 
We couple equations (|2]), ([3]), (jO]), (|9]) and flTOj) together and solve the nonlinear system with 
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a numerical method similar to Alben’s work—. We compute the free vortex sheet with an 
explicit method and solve for the unknown variables including beam position, vortex sheet 
strength, pressure and voltage using Broyden’s method^. More detailed discussion can be 
found in previous work—, including how the vorticity in the free vortex sheet is generated 
at the trailing edge using the Kutta condition, and how the Cheybyshev-Lobatto nodes are 
implemented to improve the convergence of integral-type equation. 


B. Nondimensionalization 


We nondimensionalize the governing equations by: the semi-chord b = -C, the fluid 
density pf, and the internal capacitance density Cp as the characteristic length, density and 


capacitance density, respectively. We choose the velocity scale as u = 


B 


Pfb^S 


so that 


the dimensionless inflow velocity is U = U 


Pfh^S 

B 


and f/^ represents the ratio of the fluid 


kinetic energy to the solid elastic energy. Therefore, we have the characteristic time as 


T = — = \ — -, and the following dimensionless variables: 

u \ B 


T C r t ^ ^ b^S 

C = ^, t = -, « = t = t-, p = p—, 




B 


B 


Here, / is the beam vibration frequency. We list the dimensionless parameters and the 
corresponding physical meanings in Table HI 

We then obtain the following dimensionless equations: 


MduC{s,t) = dsiTs) - ds{dsKh) + /SudsidsYfi) - \p\h 


( 11 ) 


- diV V - _ 

dttY H — z — h ^ + l3Mdiii^ = 0 
r I 


( 12 ) 
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TABLE I. Dimensionless parameters governing the piezoelectric flexible beam 


Dimensionless group 



Physical meaning 

Mass ratio 

M 

rus 

PfbS 

Ratio of the solid inertia to fluid inertia 

Reduced velocity 

U 


Ratio of the fluid kinetic energy to the 




solid elastic energy 

Piezoelectric/Mechanical 

f^M 


Ratio of the stored electrical energy in 

coupling factor 



PZT to the solid elastic energy 


Reduced resistance density 


B 


rCr, 


Pfb^S 


Ratio of the produced energy dissipated 
in resistance and stored as elastic energy 
to fluid kinetic energy 


Reduced inductance density I 


IcpB 

Pfb^S 


Ratio of the produced energy dissipated 
in inductance and stored as elastic en¬ 
ergy to fluid kinetic energy 


III. RESULTS 

A. Validation of flutter velocity and frequency without piezoelectric effect 


We first briefly review the flutter boundary for flags in an unbounded flow without the 
piezoelectric effect. The flag deformation grows exponentially when the speed of the back¬ 
ground flow exceeds the flutter speed^i. The flutter boundary curve gives the parameter 
values at which the net damping of the beam becomes zero. In the simulation, it is obtained 
by imposing small initial deflection and tracking the deflection’s growth or decay in time. 
Below the curve that represents the flutter boundary (in U), the initial perturbation on the 
beam will decay in time to zero due to the net positive damping; above the curve, small 
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perturbations will grow exponentially due to net negative damping. 


In figure [2] (a), we show the reduced fiutter speed U for beams for a large range of mass 


ratios M together with the previous experimental, numerical and analytical results, where 


all are obtained without the piezoelectric effect, 

(a) (b) 




FIG. 2. (a) Flutter velocity U at different mass ratio M. (b) Response frequency / at the flutter boundary 
with different mass ratio M. Current inviscid plate-vortex sheet method (solid line). Experimental results 
by: Watanabe et al^ conducted in air at Re 100000 (•); Huangl conducted in air at Re 40000 (■); 
Shelley et ah^ conducted in water at Re ~ 60000 (t). Numerical simulation results of viscous flow model 
by: Akcabay and YoungiS, at Re = 840 (-I-); Connell and Yuei^ at Re = 1000 (x). Numerical simulation 
results of inviscid flow model by: Tang and Poidoussisi^ with potential flow theory (dotted-diamond line); 
Michelin et aL— with point vortex model (>); Doare and Michelirii^i with double wake model (O)- Inviscid 
analytical results by: Shelley et ah^ with infinite flag model (dashed line); Eloy et al^ with double wake 
model and linear stability analysis (dash-dot line). 

In general, our inviscid vortex sheet model consistently underestimates the critical fiutter 
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speed compared to the viscous and the experimental results. This is not only due to the 
viscous effect, but due to the fact that we neglect the material damping in the model. 
Typical beams have a material damping coefficient of 1% — 5%^, and previous work shows 
that increases in material damping will lead to a delay in flutter, particularly for large M—. 
The vortex sheet model agrees well with other inviscid model results, especially in the low 
mass ratio regime [M < 1). When M is large, the beam flutters at lower mode in the 
flapping state, and the assumption that its shape is close to a travelling wave in Shelley’s 
inhnite flag modeU^ is violated. Therefore, we observe a large discrepancy between the solid 
line (current model) and the dashed line (Shelley’s model) for M > 1 in panel (a) if figure 
[21 We also notice differences in the flutter boundary between our model and the inviscid 
results from Eloy et al.^ and Tang and Poi'doussisl^. Eloy et al. used a “double wake” 
model, which has an artificial upstream wake, while our model only includes the wake shed 
from the trailing edge. This leads to different flow descriptions, especially for heavier beam 
M > 1. Tang and Pai'doussis used a vortex panel method and modeled the vortex wake as a 
series of vortices whose strength depends on the Kutta-condition. They also included 0.4% 
material damping in the model, which could be the reason that their flutter boundary curve 
is closer to the experimental results by Huang^. 

In panel (b) of figure 01 we show the flapping frequency / of the beam at the flutter 
boundary. The velocity is chosen to be right above the flutter boundary according to panel 
(a). Our results agree well with previous numerical and experimental results. 

In the flapping state, the flag shows different mode structures which are determined by 
the most unstable mode. This phenomenon is observed in experiments^*^, simulationsf^^— 
and linear stability analysis^*^. In figure [3l we show the typical mode shapes of the beam 
oscillations just above the flutter boundary. If we define the mode shape by the number 


13 


of necks in the flags’ envelops similarly to the work of Michelin et al.—, we obtain flapping 
modes ranging from 2 (one neck) and 3 (two necks) to higher values depending on the mass 
ratio M. When M is small, the relative dominance of fluid inertia and high beam flexibility 
at the flutter boundary allow higher wavenumber bending. This observation is in agreement 
with previous work-^i^il^. 


(b) 





FIG. 3. Flag snapshots corresponding to different modes of beam oscillations just above the flutter bound¬ 
ary. (a). Second mode shape, M = 4 and U = 1.4. (b). Third mode shape, M = 0.6 and U = 6. (c). 


Fourth mode shape, M = 0.2 and U = 12. 


B. Influence of piezoelectric material on flutter 

Now we consider the effect of the piezoelectric material on flutter. With the piezoelectric 
material, the effective system stiffness is increased by decreasing the fluid disturbing moment 
via the electro-mechanical coupling term with /3m in equation ([T]). This term dissipates fluid 
kinetic energy through the output circuit. Overall, the piezoelectric material functions as a 
stabilizer to the system by increasing the effective stiffness and damping. 
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1. Piezoelectric effects without inductance 


We first consider a simpler circuit without the inductor. In this case, the circuit equation 
becomes: 

diV + — + ^Mdth = 0 (13) 

r 

In hgure 01 we plot the critical flutter speed and the corresponding response frequency right 
on the flutter boundary with a hxed value of f and various values of the electric-mechanical 
coupling factor /3m over a wide range of mass ratio M. We dehne the critical flutter speed 
to be Uq and the flutter frequency to be /o when the piezoelectric material is not included, 
and plot them with solid lines in hgure 01 The critical hutter speed increases when Pm 


(a) (b) 




FIG. 4. Reduced flutter speed (a) and frequency (b) results for fixed r = 1 and /3m = 0.2, 0.5,1. The flutter 
boundary results without piezoelectric effect are plotted with solid line. 

increases for a given f. The piezoelectric material not only affects the hutter boundary, but 
also changes the vibration frequency at the hutter boundary. As the coupling factor Pm 
increases, the vibration frequency increases as shown in hgure 0](b), due to the increase in 
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effective stiffness by the piezoelectric material, which will be shown later when discussing 
the effect of the resistance. 


A flapping flag without the piezoelectric effect changes from a pre-flutter to a periodic 
flutter state to a chaotic state as the incoming flow velocity increasesl^. Similarly, with 
the piezoelectric effect, the flag undergoes periodic vibration within certain ranges of flow 
velocity. In the periodic state, we dehne the amount of the electric power generated in the 

T2' 


circuit by a time and space averaged power over one period as Pr = figure [5]^a), 

we show the power obtained right above the critical flutter speed U with a fixed value of 
f = 1 for varying M and ^m- In general, for the same mass ratio, the power increases with 
/3m due to higher electro-mechanical coupling, i.e., more fluid kinetic energy is converted to 
electric energy. The power Pr scales as as shown in figure |5]^b), where we divide the 
power by and find a collapse of all the power curves. Equation (IT^ can be rewritten as 
C = / ehm. If we assume the curvature only weakly depeuds on the piezoeleetr.e 

coupling factor, then V scales with /3m and the power scales with /3|^. As the mass ratio 
decreases, we obtain multiple peaks in Pr, and each of the peaks corresponds to a transition 
to a higher mode shape. A larger critical flutter speed is required when the flag changes to 
a higher mode shape, and the higher speed leads to more fluid kinetic energy converted to 
electrical energy, and hence larger power. The deviation in trend observed in figure E^a) for 
/3m = 1 at M = 0.2 and 0.3 is because the flag is still in the third mode shape for /3m = 1, 
but changes to a fourth mode shape for /3m = 0.2 and 0.5. Therefore a higher energy is 
obtained at low /3m values. In figure EI^c), we show a diagram of different mode shapes 
obtained with various M and /3m- For the same mass ratio, the increasing /3m increases the 
effective stiffness of the flag. Therefore, the flag oscillates in a lower mode shape, and the 
curve which differentiates different modes tends to lower mass ratios as /3m increases. 
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FIG. 5. (a). Averaged power Pr vs . mass ratio M for fixed f = 1 and various Pm = 0.2, 0.5,1- (b). 
The averaged power obtained with different Pm are rescaled by p\j and a good collapse is obtained, (c). 
Diagram of mode shapes for different mass ratio M and Pm = 0,0.5,1, 3. Pm = 0 indicates the case without 
piezoelectric effect, (d). The averaged efficiency rj vs. M for for fixed f = 1 and various Pm = 0.2,0.5,1. 

Similar to some previous work— >2^4^, we consider the energy harvesting efficiency as the 
ratio of the output electric power to available fluid power: 

{VVR) Pr 
^ PfU^AS U^A 



17 


































where A and A are the dimensional and dimensionless maximnm vibration amplitnde of 
the flag dnring one period. We show the corresponding efficiency obtained with the same 
parameters in fignreEl^d). In general, the efficiency decreases as M decreases. Althongh the 
power increases when the mode shape changes, the ffntter speed U also increases. Therefore, 
the valne of efficiency, which scales as U~^, depends on the competition of the power and 
the ffntter speed, and can either becomes larger (M = 0.4) or smaller (M = 1) when the 
mode shape changes. 

Next, we consider the effect of the resistance on the ffntter bonndary. In fignre[Hl we show 
the ffntter bonndary resnlts for a fixed valne of (3m = 0.5 and varions f. The dimensionless 
parameter f is a property of the electric circnit. When f is small, the circnit is close to a 
short circnit and the effect of piezoelectric material decreases. As shown in panel (a), the 
critical ffntter velocities are closer to t/o as f —)■ 0. When f is large, the circnit is close to an 
open circnit. In the limiting case when f —)■ cxo, the circnit eqnation can be rewritten as 

V = -I3mk (15) 

therefore, the flag behaves eqnivalently to a flexible flag with a higher bending rigidity of 
i?(l + /S£f) and the critical velocity shonld be Uq (l + (3l^) . The piezoelectric material 

therefore acts to increase the effective bending stiffness of the system. When the resistance 
has a moderate valne, the voltage does not depend linearly on the cnrvatnre. However, it 
can still be represented by the prodnct of (3m and a fnnction of cnrvatnre, which snggests 
an increase in the system’s bending stiffness with an increasing (3m- We notice that this 
observation is consistent with the resnlts obtained by Doare and Michelin^. In fignre [Hl(b), 
we plot the stability bonndary resnlts for fixed (3m = 0.5 and f = 1,10,100 together with 
the curve Uq (l + and Uq, and we can see that when f = 100, the curve approaches 

the curve t/o (l + (3\j') . The corresponding vibration frequencies for varying f are shown 
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(b) 




(c) (d) 




FIG. 6. (a). Critical flutter speed for fixed /3m = 0.5 and r = 0.01,0.1,1- The flutter speed without 
piezoelectric effect is Uq plotted with solid line. (b). Critical flutter speed for fixed Pm = 0.5 and r = 
1,10,100. Solid line shows Uq and dash-dot line shows the curve Uq (l -I- /3^) . (c). Frequency values at 

the flutter boundary for fixed Pm =0.5 and f = 0.01, 0.1,1. The frequency without piezoelectric effect /g is 


plotted with solid line. (d). Frequency values at the flutter boundary for fixed Pm = 0.5 and r = 1,10,100. 

_ 1 / 

Solid line shows /g and dashed-dot line shows the curve /g (l + Pm) 
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in panels (c) and (d) of figure 01 Similarly to the critical speed, when r —)■ 0, the frequency 

_ -| j r\ 

tends to /o; when f —)■ cx), the frequency tends to /o (l + 

Unlike the trend with varying /3 m, the flutter velocity as well as the flutter frequency 
change non-monotonically with f, as shown in hgure [HI In hgure [71 we plot the flutter 
velocity, the corresponding vibration frequency, the averaged electric power and efficiency 
right above the ffutter boundary varying over a range of f with a hxed mass ratio M = 10 
and /3 m = 0.2, 0.5,1. We also plot Uq and /o with solid lines for comparison. 

For all three Pm values, the nondimensional resistance which leads to largest flutter speed 
is obtained around f = 1, when the fluid kinetic energy is balanced by the energy dissipated 
in the resistor and stored as electric energy in the piezoelectric material. The maximum 
flutter frequency and averaged electric power are also achieved at these intermediate points. 
Since the efficiency scales with the inflow velocity as U~^, the largest rj is obtained at a 
smaller value of f as shown in panel (d). The exitance of an optimal resistance load was 
studied previously for piezoelectric beams in a vacuum with base excitations22>^, and was 
also observed for piezoelectric beams in viscous flowl^. 

Another important observation from hgures 0| and [H] is that the numerical results show 
that the increase in flutter speed by PZT (with increase in Pm and f) is more rapid for 
cases with higher M [M > 1), i.e., heavier flags in a lighter fluid. For lighter flags in a 
heavier fluid (M -C 1), the piezoelectric material will have a negligible impact on the beam 
vibrational response for cases with /3 m < 0.2 or f < 0.1. 

2. Piezoelectric effect with inductance 

Now we consider the effects of inductance on the flutter boundary. Similarly to the effect 
of the resistance, when inductance tends to zero (/ —)■ 0), the circuit tends to a short circuit 
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(a) (b) 



FIG. 7. (a). Flutter velocity U vs. resistance f. The solid line without symbols denotes Uq, the flutter 
velocity when no piezoelectric material is included, (b). Vibration frequency right above the flutter boundary 
/ vs. resistance f. The solid line without symbols denotes /o, the frequency when no piezoelectric material 
is included, (c). Averaged power Pr vs. r. (d). Averaged efficiency rj vs. f. Other parameters used for 
three panels are M = 10, Pm = 0.2, 0.5,1. 

and the effect of the piezoelectric material is negligible. When indnctance tends to infinity 
(/ —)■ cxd), the circnit tends to an open circnit and the circnit is a pnrely resistive circnit as 
discnssed in the previons section. 
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When the inductance has a moderate value, the parallel RL circuit has a natural frequency 
of fe = cases, we find that the piezoelectric flag always exhibits a chaotic 

motion in the post-flutter regime. Periodic solutions with small perturbations are found 
close to the flutter boundary when < 1- Therefore, instead of considering the power and 
efficiency over one period, we use an averaged power Pr and efficiency rj over a long time in 
the large-amplitude regime as measurements. In hgure El we show the critical flutter speed, 
vibration frequency and the corresponding averaged power and efficiency over a range of the 
circuit frequency fe with hxed mass ratio M = 1, f = 1 and various /3m = 0.2, 0.5, 0.8. 

In panel (a) of hgure El we notice a drop in the critical hutter speed as the circuit 
frequency fe increases above a certain value, and the hutter speed after the drop is smaller 
as /3m increases. The hag is able to hutter at a much lower speed due to the “lock-in” 
between the circuit resonant frequency and the hag resonant frequency, i.e. / = fe- In panel 
(b), we plot the vibration frequency / and the circuit frequency fe together with a solid 
line to denote the resonance situation where f = fe- We note that since the hag undergoes 
chaotic motion in this state, the oscillation frequency of the hag is not single-valued. We 
instead use the peak of the power spectrum of the frequency in the large-amplitude happing 
state. Therefore, we observe some discrepancies between the hag vibration frequency and 
the circuit resonant frequency. Nonetheless, in a certain range of I or fe, the slope of the 
two frequency curves match well with each other, and the hutter speed decreases rapidly 
in the “lock-in” range. This “lock-in” phenomenon was also observed by Xia et alM, and 
they also noticed a sudden decrease in the hutter speed with a slender body model. In the 
“lock-in” regime, we notice a rapid increase in the power dissipated in the resistor, as well 
as the efficiency, as shown in panels (c) and (d). This is an important improvement for an 
energy harvester, since one of the limitations of energy harvesters is their requirement for 
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FIG. 8. Flutter boundary results for M = 1, r = 1, and /3 m = 0.2,0.5,0.8. (a). Critical flutter speed U vs. 
circuit frequency fe = — =rF}- Tbe solid line indicates the flutter speed Uq when no piezoelectric material 
is included, (b). Vibration frequency / vs. circuit frequency fe- The solid line shows the line when / = fe, 
which corresponds to the “lock-in” situation, (c). Averaged power dissipated on the resistance Pr vs. fe- 
(d). Averaged efficiency rj vs. fe- 

flow speed which is difficult to achieve in a natural environment. By tuning the inductance 
into the “lock-in” regime, the energy harvester is able to flutter at a much lower speed and 
produce a much higher energy. 
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The flutter speed is greatly reduced in the “lock-in” regime as the flag jumps to a lower 
mode shape which requires less fluid kinetic energy to initiate flutter. If the flag is already 
in the lowest mode shape (second mode shape) and cannot be reduced further (such as for 
M > 1 according to figure [21(a)), we And that the flutter speed will not be reduced even 
though flag is in the “lock-in” condition. In figure [9l we show the flutter speed and the 
corresponding vibration frequency over a range of circuit resonant frequencies with fixed 
values of /3 m and f and different mass ratios. In figure (TUI we show the snapshots of the flag 
at different mass ratios with and without the inductor. The values of I are chosen so that 
the vibration frequency of the beam is “lock-in” to the circuit resonant frequency. 

(a) (b) 




FIG. 9. (a) Flutter speed U vs. fe = —(b) vibration frequency / vs. /e at M = 0.1,1,10, /3 m = 0.5, 
and f = 1. Dashed lines denote the flutter speeds and frequency without piezoelectric effect for different 
mass ratios. Solid line in panel (b) denotes the curve / = fe- 

Without the inductor, the flag flutters in the 5th and 3rd modes respectively for M = 0.1 
and 1, as shown in figures [TOf d) and (e). However, in the “lock-in” regime, the flag flutters 
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FIG. 10. Flag snapshots at different M and fixed values of /3 m = 0.5 and f = 1, with the inductor (a)-(c) 
at “lock-in” / = /e; and without the inductor (d)-(f). (a). Second mode shape, M = 0.1, U = 5.8, and 
I = 0.04. (b). Second mode shape, M = 1, C7 = 4.2, and I = 0.02. (c). Second mode shape, M = 10, 
tJ = 1.9 and I = 0.6. (d). Fifth mode shape, M = 0.1 and U = 19.5. (e) Third mode shape, M = 1 and 
ty = 5.6. (f). Second mode shape, M = 10 and C/ = 1.9. 


in the 2nd mode for both cases as shown in panel (a) and (b). When /e is small (J ^ 0), 
the circuit is close to a purely resistive circuit, and the flutter speed and vibration frequency 
with the PZT are larger than Uq and /q. Notice in hgure|l](a), the piezoelectric material is 
more effective for larger M [M > 1) at /3m = 0.5. Therefore, the flutter velocity U increases 
more for M = 10 than the other two cases. As I decreases or j/ increases and approaches 
the lower bound of “lock-in” regime, the flutter speed and vibration frequency of the flag 
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drop rapidly at first; as fe continues increasing, the flutter speed and vibration frquency 
gradually increases, eventually back to Uq and /o, as T ^ 0 or /e —)■ oo corresponds to the 
short circuit scenario. When M = 10, although the two frequencies still match each other 
around /e = 0.2 (J ^ 0.6), the flutter speed simply decreases to Uq as /e increases. In 
panel (c) and (f) for M = 10, the flag flutters in the 2nd mode with/without the inductor, 
although the motions are different. In addition, notice in £gure[9]that the “lock-in” range 
is wider for lower M, where there is more room for the vibrating flag to drop down to the 
lowest flutter mode (2nd mode). This indicates that when the fluid is light or the solid is 
heavy (M > 1), the RL shunting circuit is not effective at modifying the flutter velocity or 
frequency. On the other hand, the results indicate that the RL shunting circuit can be used 
to tune the vibration responses of a lightweight piezoelectric flag in a heavy fluid such as 
water (M < 1). 

IV. SCALABILITY OF PIEZOELECTRIC BEAMS 

Now, we consider several commonly used commercial materials and show how to apply 
the above discussion to actual devices in different fluids. We hrst show some typical values 
of the above dimensionless parameters for PZT plates in air and water. We consider a 
rectangular plate with chord length of 0.05 m, a thickness-to-chord ratio of 0.5%, and a 
span-to-chord ratio of 1. These are typical dimensions for a thin plate and are similar 
some of the experimental systems^'ii'^. PZT-5II is selected to be the piezoelectric material 
due to its high piezoelectric coupling coefficient relative to other PZT material. Each PZT 
patch is assumed to have a hxed thickness of 0.02 mm, and the substrate material has 
a fixed thickness Qs of 0.06 mm. The dimensional piezoelectric coupling coefficient is then 
6.6 X 10“®m- N/V, and the internal capacitance is 6.4 x 10“® F for each patch. The loaded 
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resistance is chosen to be 1 kf2 and the inductance is 0.1 H. 


In Table |TT1 we list the density, Young’s modulus and Poisson’s ratio for four common 
substrate materials (Balsa Wood, Polymer (POM), Aluminum Alloy and Steel) spanning a 
wide range of material density and elasticity values. We also list the material properties of 
the PZT-5H. 


TABLE II. Parameters of the Solid Material and PZT-5H 



Balsa 

Wood 

POM 

Alumina 

Alloy 

Steel 

PZT-5H 

Density p (kg/m^) 

97.1 

1480 

2800 

7850 

7700 

Young’s Modulus E (GPa) 

1.8 

3 

79 

210 

66 

Poisson’s Ratio 

0.36 

0.35 

0.33 

0.3 

0.34 


The bending rigidity for the composite piezoelectric flag is calculated by: 


B = El = 


Es qlS , E, 


+ 


;S{-qlqp + qsql + ^ql) 




(16) 


where Eg, Vg and Ep, Vp are the Young’s modulus and Poisson’s ratio for the substrate and 
the PZT patch, respectively. In hgure [TTl we show the dimensionless M, Pm and f for 
different materials according to the dehnitions given in Table [H We consider two types of 
fluids: air and water to cover the range of M, i.e., heavy flags in light fluid and light flags 
in heavy fluid, respectively. 

First, we note that (3m is unchanged as the fluid medium changes from air to water since 
it is not a function of p/, as shown in Table [H Notice from Table m the PZT patch has a 
large density and Young’s modulus compared to many substrate materials. In general, the 
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FIG. 11. Range of dimensionless parameters corresponding to different materials. The open symbols 
indicate materials in air, and the filled symbols indicate materials in water. o/v/A/D: steel/Aluminum 
alloy/POM/Balsa wood in air; •/t/a/B: steel/Aluminum alloy/POM/Balsa wood in water, (a). M vs. 
Pm', (b). M vs. f; (c). M vs 1. 

mass ratio of various flags in water is much smaller than those in air, M is around 10 — 30 
for flags in air, and therefore they mainly flntter in the second mode. For flags in water, M 
is aronnd 0.01 — 0.03; they flntter in a higher mode and the critical flntter speed withont 
PZT is mnch larger, as shown in hgnre [2l 

The dimensionless resistance density f and indnctance density I both decrease snbstan- 
tially when the flnid becomes heavier as shown in panels (b) and (c) of hgnre [HI We 
also note that the dimensionless f and I depend linearly on the valnes of resistance and 
indnctance, and therefore can be easily adjnsted by changing their valnes. 

Now, we take POM as an example and show the dimensional hntter speed, vibration 
freqnency, averaged power and efficiency in air and water nnder different conditions in Table 
m We consider both the resistance only circnit and the RL circnit. The effectiveness of 
the resistance only circnit on modifying the hntter speed and freqnency is larger for hags in 
air, and practically negligible for hags in water. By nsing the optimal resistance valne in a 
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resistance only circuit shown in Table IIIIl the flutter speed is increased by 50% compared 
to the case without PZT when the flag is placed in air, but only by 2% when it is placed 
in water. A larger electric power and efficiency are also obtained with the optimal R when 
the flag is placed in air. For heavy fluids, the optimal R value which leads to a larger flutter 
velocity results in a smaller power and efficiency. The voltages for different R values are 
almost identical in this case, so the power and efficiency decrease with a higher resistance. 

On the other hand, for an RL shunting circuit, the flag in the lighter fluid mainly flutters 
in the second mode (the lowest flutter mode), and hence the flutter speed and the frequency 
will not be greatly reduced by the inductance, and there is not an optimal inductance value. 
The RL circuit is more effective for lighter flags in a heavier fluid, as shown in the example 
of water in Table m We consider the RL circuit with the optimal R value. When the flag is 
placed in water, by tuning the inductor to the “lock-in” range, the flutter speed is reduced 
by 75% compared to the case without PZT. The averaged power increased by a factor of 
10 compared to the case with a purely resistive circuit and the efficiency increased by a 
factor of 1000. Note that both circuits have important practical applications. An increase 
of the flutter speed is typically desired for vibration control applications, while an decrease 
is typically desired for energy harvesting applications. 


For an actual device with a larger size and complete geometric similarities, i.e., if the 
thickness-to-chord ratio, thickness ratio of the PZT patch to substrate material, and span- 
to-chord ratio are fixed, then the dimensionless M, /Sm and f are unchanged with the length 
scale for the same PZT and substrate material, fluid and resistance values. According to 
Table [H /3m = /3m \ -• If we assume the length scale is b, then the dimensional /3m 

V 

scales as b^, B scales as b^, and capacitance density is unchanged with the length, which 
implies that /3m is unchanged with the length scale. A similar analysis shows that f is also 
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TABLE III. The effect of the PZT on the flutter speed, frequency, power and efficiency for a POM 


flag in air and water 


Fluid 

Air (pf = 1.18 kg/m^) 

M = 13.6, /3m = 0.26 

Water 

M = 

{pf = 998 kg/m^) 

0.016, /3m = 0.26 

Resistance only 

No PZT 

R 

Optimal R 

No PZT 

R 

Optimal R 

circuit 


lOOOH 

127H 


lOOOH 

3477H 

Flutter speed (m/s) 

14.9 

20 

22.4 

18 

18.2 

18.4 

Frequency (Hz) 

131 

137 

118 

323 

150 

150 

Power Pfi (W) 

- 

0.67 

3.55 

- 

0.93 

0.27 

Efficiency r] 

- 

0.023 

0.078 

- 

3.4e-4 

l.Oe-4 

RL circuit 

No PZT 

L 

No Optimal 

No PZT 

L 

Optimal L 


- 

O.IH 

L 


O.IH 

0.93H 

Flutter speed (m/s) 

14.9 

14.9 

- 

18 

8.7 

4.9 

Frequency (Hz) 

131 

131 

- 

323 

44 

13 

Power Pr (W) 

- 

0.32 

- 

- 

1.78 

3.1 

Efficiency r] 

- 

0.0015 

- 

- 

0.02 

0.14 


invariant with the length. I is proportional to 6, but when the length increases, the same 
dimensionless value can be obtained by changing the value of the inductor. Therefore, the 
dimensionless results obtained from the model problems can be directly applied to the full 
scale problems. The dimensional velocity is invariant with respect to b, the frequency scales 
as b~^ and the power scales as 6^. Therefore, when the problem is extended to a larger size, 
the flutter speed is unchanged, but a smaller frequency and a larger power are expected. 

On the other hand, the actual thickness of the PZT patch can not exceed several cen- 
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timeters due to material properties and manufacture restrictions. Thus, the thickness ratio 
of the PZT patch and the substrate material can be much smaller than that of the model 
scale problem, which leads to a smaller [3m- Another important issue is that the thickness to 
chord ratio is generally larger for an actual device, otherwise the plate may undergo material 
failure, which is not considered in the current analysis. The bending rigidity B is larger for 
a thicker plate, and it leads to a larger flutter speed which scales as according to Table 
m Moreover, although the substrate material can be scaled, the material properties of the 
PZT patch can be difficult to scale. For example, very large PZT patches are generally 
not feasible in reality, so several discrete patches are often linked together as an alternative. 
These issues along with the stability of the particular material in a certain fluid must be 
considered when extending the model scale results to full scale problems. 

V. CONCLUSIONS 

In this work, we investigate the effect of the PZT material on the flutter speed, vibration 
mode and frequency, and energy harvesting power and efficiency of a flexible flag in differ¬ 
ent fluid medium. We develop a fully coupled ffuid-solid-electric model by combining the 
inviscid vortex sheet model with a linear electro-mechanical coupling model. A cantilever 
flag with two layers of PZT and a uniform inflow is studied as a model problem. Based on 
a nondimensional analysis, four dimensionless parameters M,(3m,t,1 were found to govern 
the reduced flutter velocity U and frequency / of the piezoelectric flag. 

When the mass ratio M is small, the flag flutters at higher mode, as the relative high fluid 
inertia and fluid disturbing force allow higher wave number deforming patterns. Flags with 
lower M also require a larger inflow velocity to initiate flutter, due to the higher relative 
fluid inertia and the damping resistance against motion. 
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Two types of output circuits are considered in this work: a purely resistive circuit and a 
resonant RL circuit. For a purely resistive shunting circuit, the PZT material functions as a 
stabilizer to the system by increasing the effective system stiffness and damping. Therefore, 
as /3m increases, the flutter speed, as well as the vibration frequency, averaged power and 
efficiency also increase. In particular, we hnd that the power scales as For a resistance 
only circuit, the optimal resistance to maximize the flutter speed and the power is obtained 
at near f = 1, when the fluid kinetic energy is balanced by the energy dissipated in the 
resistor and stored as electrical energy in the piezoelectric material. The results also show 
that a purely resistive circuit have a greater impact for heavier flags in lighter fluids, i.e., 
M > 1. For M 1, the PZT flag with a resistance only circuit has a negligible impact if 
Pm < 0.2 and/or f < 0.1, which is the case for most PZT flags in water. For a resonant RL 
circuit, an optimal inductor is found when the circuit frequency matches the flag’s vibration 
frequency. In this “lock-in” regime, the flutter speed and frequency are greatly reduced and 
the flag will flutter in a lower mode. The “lock-in” range is found to be larger, and the 
flutter speed and the vibration frequency drop more in this range for lower M. Maxima 
electric power and efficiency are also obtained in the “lock-in” regime. 

We also consider a model scale set-up with several commonly used commercial material 
in air and water. Typical ranges of dimensionless parameters are obtained for four types of 
materials that spans a wide range of solid density and rigidity values. For POM, we hnd 
that the resistance only circuit is more effective when the hag is placed in a lighter huid like 
air. An RL circuit is able to signihcantly reduce the hutter speed and frequency when the 
hag is placed in a heavier huid like water, but is not as ehective for similar hags in air. For 
vibration control applications, an increase in the hutter speed is typically desired. On the 
other hand, for energy harvesting applications, a lower hutter speed is desired. 
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An extension of this work is to consider the application of the PZT material on passive 
vibration control and for energy harvesting with a vibrating foil, which requires generalizing 
the current vortex sheet model to more complicated geometric shapes, and considering a 
combined bending and twisting deformation. Moreover, the influence of viscous and 3-D 
effects on electro-fluid-solid interactions requires additional study. 
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